On the motifs distribution in random hierarchical networks 
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The distribution of motifs in random hierarchical networks defined by nonsymmetric random 
block-hierarchical adjacency matrices, is constructed for the first time. According to the classi- 
fication of U. Alon et al of network superfamilies [llj by their motifs distributions, our artificial 
directed random hierarchical networks falls into the superfamily of natural networks to which the 
class of neuron networks belongs. This is the first example of "handmade" networks with the motifs 
distribution as in a special class of natural networks of essential biological importance. 



I. INTRODUCTION AND BASIC DEFINITIONS 



Most commonly, the hierarchy of states emerges in 
many-particle systems of various origins with a large 
number of "frozen" constraints with different scales, 
which generate multidimensional hypcrsurfaces of poten- 
tial energy (or free energy) with an astronomically large 
number of local minima. Typical examples of such sys- 
tems (often referred to as complex systems) are glasses 
and globular proteins. The hierarchical concept applied 
to such systems presumes that local minima of the en- 
ergy landscape are clustered into hierarchically embed- 
ded basins of minima. Namely, each large basin consists 
of smaller basins each of which in turn contains embed- 
ded still smaller basins, and so on. Local minima basins 
are separated from one another by hierarchically ordered 
barriers (the smaller the basins, the "lower" the barriers 
separating them). 

Apart from the dynamic contents of the hierarchical 
concept, the determination of the hierarchical organiza- 
tion of "ultrametric phase spaces" in the observed sta- 
tistical regularities is of considerable interest. A visual 
example of the such a structural organization is the so- 
called crumpled globule discussed for the first time in 
[l|. The thcrmodynamically equilibrium spatial confi 
uration of such a globule resembles the Peano curve 
embedded into a 3D space. Spatial packing of a crum- 
pled globule can be represented schematically by a single 
folded motive reproduced on a growing scale. The hi- 
erarchical packing naturally leads to a block-hierarchical 
network of contacts between the links of a chain described 
by a block-hierarchical matrix of contacts. Naturally, the 
presence of inhomogeneities in the hierarchy of crumples 
introduces randomness in the block-hierarchical network 
of contacts, which requires the determination of statis- 
tical characteristics of an ensemble of random block- 
hierarchical matrices of contacts. In recent works [3, |4j 
we have considered statistical properties of random hier- 
archical networks defined by adjacency matrices in form 
of block-hierarchical Parisi matrix [5|. Remind that a 



network is a set of vertices (or nodes) and connections 
between them (links or edges). We suppose that in the 
network there are no any self-connections and multiple 
edges. The network is random if any link occurs with a 
certain probability. The network is directed if any link ei- 
ther has an orientation (i — > j) or is bidirectional (i •(-> j). 
Otherwise the network is non-directed. 

The investigation of statistical properties of random 
graphs and networks implies studying of spectral prop- 
erties of their adjacency matrices (e.g. 0,01), the same 
question can beenposed for block-hierarchical networks. 
It was found in [3, [I| that the spectral density of adja- 
cency matrices has power law ("heavy") tails, typical for 
scale-free networks. This observation has been supple- 
mented by direct investigations of such typical statisti- 
cal properties of networks as vertex degree distribution, 
which turned out to be abnormally wide (but not scale- 
free). Hence following the conventional classification (e.g. 
[8j), random hierarchical networks could be attributed to 
the class of scale-free (according to the spectral density) 
or polyscale (according to the vertex degree distribution) 
networks. 

Scale free- networks are associated with a variety of 
structures and systems, such as protein folding and 
biopolymer dynamics; cell metabolism; neural, informa- 
tion and communication networks; various evolutional, 
ecological, social and economical systems. Statistical 
characteristics of many natural networks are described 
in the review [8| . Because of wide usage of a "network 
paradigm" it seems quite natural that a "handmade" de- 
sign of artificial networks with some observed statistical 
characteristics close to that of natural networks is of pri- 
mary importance. Such a design might be very useful 
tool for searching for the correlations between the net- 
work structural organization and the functions. 

Until recently building of scale-free networks was 
based in almost all works on a step-by-step growing pro- 
cess based on the preferential attachment method [9( and 
its various modifications. In these approaches the new 
vertices are connected to the existing ones with proba- 



bility which depend on their current vertex degree. Most 
of the statistical characteristics of artificial scale-free net- 
works including spectral distribution of adjacency matrix 
were obtained for networks built in this way. It should be 
noted that the preferential attachment process, which re- 
alizes locally inhomogeneous vertex grouping, implicitly 
implies some mechanisms controlling the current state of 
a network with long-term "evolutionary memory" . 

Unlike the preferential attachment-like methods, 
building of hierarchical networks is based on constructing 
of hierarchically embedded clusters of links. The config- 
uration of links is usually described by an adjacency ma- 
trix A in which for matrix elements a,y one has atj = 1 
if there is a link connecting nodes i and j and a%j = 
otherwise. Adjacency matrix A for non-directed graph 
is symmetric, i.e. a^ 
directed network a 
link i — > j and a 



= 1 and 



To the contrary, for the 
for the oriented 



v ' 



ij 



'./' 



1 for the bidirectional one 
i •O- j. Therefore the generic adjacency matrix A is not 
necessarily symmetric. 

It is known that the classification of many natural 
networks according to their vertex degree distribution, 
or clustering coefficient is too rough and does not pro- 
vide any relevant information about the internal network 
structure. Much more detailed information about the 
network structure can be provided by investigating the 
local topological characteristics, the so-called motifs and 
their distributions [10|,[ll|]. For example, it is known that 
all networks, according to their three-vertex oriented mo- 
tifs distribution, can be divided into four super families 

B 

In this letter we announce the results of the investi- 
gation of motifs distribution in random hierarchical net- 
works. The key outcome consists in the fact that the mo- 
tifs distribution of random block-hierarchical networks 
clearly falls into one of the universal supcrfamilics, which 
includes, in particular, networks of neurons. Besides, we 
claim the existence of a phase transition (in respect to 
motifs' distribution) in an ensemble of block-hierarchical 
networks in the thermodynamic limit. 



II. DISTRIBUTION OF MOTIFS IN 
HIERARCHICAL NETWORKS 



Remind that local topological properties of networks, 
both directed and non-directed, for given number of ver- 
tices and vertex degree distribution can be characterized 
by the rates of connected subgraphs. Since the number 
of such subgraphs grows combinatorially with their size, 
usually only small subgraphs are considered. In partic- 
ular, in the works [ljj, [llj only subgraphs of size 3 (tri- 
ads) were analyzed for directed networks. There are 13 
different configurations of such triads. They all are enu- 
merated in the Fig[TJ 



1 2 3 4 5 6 7 8 9 10 11 12 13 
Figure 1: Connected subgraphs-triads for directed networks. 



The rates of subgraphs in a given network depend on 
the vertex degree distribution. This complicates the com- 
parison of networks of different sizes and different vertex 
degree distributions by the rates of their subgraphs. In 
order to compensate these differences, the procedure of 
so-called network randomization was proposed in works 
jlClL lll| . In this procedure the network experiences multi- 
ple permutations of links under the condition of conserva- 
tion in each vertex of the number of incoming, outcoming 
and bidirectional links. Using this method an ensemble 
of randomized versions of a given network is generated, 
and for every subgraph the statistical significance 



Z k = 



N k - (N k ) 



rand 



<Jk 



(1) 



is calculated, where N k is the amount of fc-th subgraphs 
in the initial network and (Nk) rand and o~ k are corre- 
spondingly the mean and the standard deviation of N k 
for the randomized networks. Subgraphs with the statis- 
tical significance essentially exceeding 1 are called motifs 
jlfj |. The motifs' distribution of the network under con- 
sideration is characterized by a significance profile which 
is a normalized vector 



P = {Pl,-,Pm} 



(2) 



of statistical significance for all subgraphs of given size. 
The components of the vector p are: 



Z k 



Pk = 




(k = l,...,m) 



(3) 



It has been demonstrated in the papers [lCj, [ll| that 
significance profile distribution could be used to divide 
networks into superfamilies. For directed networks only 
4 of such superfamilies were determined. The networks 
with considerably different functional properties, for ex- 
ample, the neuron networks and transcriptional networks 
in unicellular organisms, belong to different superfami- 
lies. In [10|, LLL[ also the undirected networks were clas- 
sified according to their tetradic motifs. These networks 
were separated in four superfamilies as well. 

It is interesting that artificial random scale-free net- 
works generated by the preferential attachment method 
form a separate supcrfamily which does not coincide with 
any superfamily of real networks. In the Figf2] we have 
reproduced from [1JJ the significance profiles of one of 
the superfamilies for directed networks, which later on 
will be compared with our results on the distribution of 
motifs in the block-hierarchical random networks. 



Experimental data (reproduced from U.AIon et al. [1 1]) 
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Figure 2: (Color online) Significance profile distribution of 
motifs-triads corresponding to the superfamily of directed 
networks, to which the network of synaptic contacts between 
the neurons in C.elegans belong (reproduced from [lUp. 



The generic procedure of the random block hierarchi- 
cal (RBH) network construction is as follows. Taking N 
points as potential vertices of our forthcoming network, 
we raise a hierarchical network by connecting the vertices 
by edges in a specific way. We consider the adjacency ma- 
trix in form of a p-adic translation-noninvariant Parisi 
matrix A. This matrix is shown in Fig|3]for p = 2. 
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Figure 3: Random p-adic (p — 2) block-hierarchical adja- 
cency matrix. 

Since we are aimed to describe directed networks, the 
matrix oi A may not be symmetric: ay 7^ aj%. All matrix 



elements, a.; 



» 



and a-j 



of , are the Bernoulli 



{qi,q2, ■■-,qr}- Thus, the elements of A, being the ran- 
dom variables, are hierarchically organized in probabili- 
ties. In case of directed networks the matrix elements 
above and below diagonal were generated independently. 
Below we consider the set of probabilities, {Q}, with 



Q~i 



= p-^< 



(/x>0) 



(5) 



In general p > 2 (we consider the case p = 2), 7 = 
l,2,...,7 max is the hierarchy level, and \i > is a pa- 
rameter. This methods allows to rise scale-free networks 
with two important features. First, formation of clus- 
ters of links on each hierarchy level 7 is uncorrelated. 
Second, random subgraphs associated with different hi- 
erarchy levels could be different, so the whole network is 
not necessarily homogeneous. 

The systematic study of statistical properties of ensem- 
bles of random graphs (networks) deals with the inves- 
tigation of the spectral properties of a graph adjacency 
matrix jg,0]- Let A^ (1 < i < N) be the eigenvalue of the 
adjacency matrix. The spectral density of the ensemble 
of random symmetric adjacency matrices is defined in the 
standard way, 



1 N 



{qi,Q2,---,gr} 



where 



denotes the 



(0) 



the 



) {qi.q- 2 ,...q n } "cuuucs une averaging c 
distributions of the matrix elements of A. Comput- 
ing numerically the spectral density, p(A), of networks 
with symmetric block-hierarchical adjacency matrices, 
we found that the tails of the spectral density p{\) follow 
a power-law ( "heavy tail" ) asymptotic behavior 

p(\) ~ |A|~* (7) 

In Q we have found that for fi e]0, 1[ the exponents x(^) 
takes the values slightly below x = 2- From this point of 
view the random hierarchical networks with symmetric 
adjacency matrices are scale-free. 

The hierarchical structure of clusters leads to a scale- 
free (polyscale) structure of network in a broad range of 
the parameter \x (see 0] ) . This fact suggests that random 
hierarchical networks (both symmetric and nonsymmet- 
ric) might serve as a model for certain families of natural 
scale-free (polyscale) networks formed without a specific 
growing mechanism. 



distributed random variables: 



{<>,&M} 



1 with the probability g 7 
with the probability 1 — q 7 



(4) 



where 7 counts the hierarchy levels (1 < 7 < 7 max = T) 
and n enumerates different blocks corresponding to a 
given hierarchy level 7 (see Fig|3|). Note that the prob- 
ability q 7 does not depend on n. The full ensemble of 
N x N matrices A, where N — p r (p ~ 2), is com- 
pletely determined by the set of probabilities, {Q} = 



III. RESULTS 



Distribution of motifs— triads in directed 
hierarchical networks 



The significance profile we have computed by the 
method described above, using triads for directed net- 
works. Since the networks are generated at random, for 



each value of fi the significance profiles are averaged over 
an ensemble of corresponding random hierarchical net- 
works. The typical distributions of motifs for directed 
hierarchical networks is shown in FigJU As one can see, 
the distributions for different values of fj, look very simi- 
lar and are not too sensitive to the concrete value of [i, 
meaning that the hierarchical networks for different /z's 
are topologically similar as well. 
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Figure 4: a) (Color online) Normalized motifs distribution for 
hierarchical 2-adic random network of size 32 x 32 for different 
fi (see the text for details); b) Normalized motifs distribution 
for hierarchical 2-adic random network of different sizes for 
single fixed value /i = 0.8. 



Comparing our distribution of motifs for directed hier- 
archical networks depicted in FigUwith the "second su- 
pcrfamily" in the classification of U. Alon et al shown in 
FigO (look [Tlj for more details) , one sees that in a broad 
range of /i's our directed hierarchical networks clearly fall 
into the "second superfamily" , to which, for example, the 
neuron networks belong. 

We would like to emphasize that the hierarchical 
directed networks is, apparently, the first example of 
"hand-made" artificial networks topologically similar to 
a certain superfamily of natural networks in terms of lo- 
cal topological properties. We should stress that hier- 
archical random networks can be built by uncorrelatcd 
generation of clusters of links, unlike the essentially cor- 
related preferential attachment procedure. In the light of 
results obtained, this feature looks particularly interest- 
ing in the context of modelling of biological operational 
systems and their evolutionary prototypes (e.g. |12|). 
From this point of view the hierarchical networks could 
be of particular interest for neuron networks modelling. 



B. Phase transition 

We have checked our distribution of motifs (Fig0|) of 
directed hierarchical network on the stability. To study 
this question, the following numerical experiment has 
been performed. First, we have generated the block- 
hierarchical adjacency matrix of some directed graph, 
and then we randomly spoiled this block-hierarchical 
structure by the following procedure. 

To be precise, we scanned once all matrix elements of 
the adjacency matrix row-by-row from the first element, 
On, t° the last one, clnn- Each element a%j we have inde- 
pendently switched to the opposite value with the prob- 
ability /, i.e. if dij = 1, then with the probability / the 
element atj can take the value and vis versa. Obviously, 
for the noise / = 1/2 after one run over all matrix ele- 
ments, we have destroyed all hierarchical blocks and have 
generated a completely random adjacency matrix corre- 
sponding to the Erdos-Renyi random graph [13[ with an 
appropriate distribution of motifs. Since our motifs' dis- 
tribution is measured off the distribution of random un- 
corrected graphs, the corresponding significance profile 
shown in Fig|4]would be definitely for / up = 1/2 for all 
13 configurations of triads. 

To characterize quantitatively the degree of similarity 
between motif's distributions for different /, we define 
the scalar product, rj{f): 



13 



V(f) = P(/)P(0) = X>(/)1*(0) 



(8) 



fc=i 



where p(/) is the distribution of motifs for a given value / 
of the noise, and p(0) is the reference (initial) distribution 
of motifs for / = 0. The value r\ serves as an "order 
parameter" and its critical behavior is shown in the Figj5] 
for three different sizes of networks 32 x 32, 64 x 64 and 
128 x 128. 
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Figure 5: (Color online) Signature of the phase transition in 
distribution of motifs. 



Varying the "intensity of noise", /, from / = up to 
/ up = 1/2, we have noticed that the order parameter r\ 



becomes statistically indistinguishable from for values 
of / significantly less than / up = 1/2. Namely assuming 
the critical behavior 



V(f)= C|/cr-/|' 



(9) 



we have found the following numerical values c w 1.72; 
f CI w 0.3; /? ps 0.43 for the best fit of the data corre- 
sponding the matrix of size 128 x 128. We can interpret 
this behavior as a signature of a possible phase transition 
characterized by the behavior of the order parameter rj(f) 
in the block-hierarchical networks in the thermodynamic 
limit. 

In order to demonstrate that the existence of the phase 
transition at finite temperature is very natural for the 
hierarchical system, we consider in the Appendix the toy 
model of the spin system on a complete graph with block- 
hierarchical coupling known in the literature as a "Dyson 
hierarchical model" [1J, [Tj|. The method used here is 
based on the so-called p-adic Fourier transformation (see, 
for example |16[ for details) and allows to re-derive the 
classical results on Dyson model in few lines. 

The physical meaning of the hierarchy of phase transi- 
tion in the hierarchical model discussed in the Appendix 
lAlis very clear. When the temperature decreases, at first 
critical temperature, T cr , the spins become correlated 
within the smallest clusters only, where the interaction is 
the most strong. However, the clusters of 2nd and higher 
levels of hierarchy remain still uncorrelatcd because the 
interaction of spins inside them is weaker. At the second 
critical temperature, T„ ' (T cr < T cr ) the smallest clus- 
ters remain correlated, together with the clusters of the 
2nd hierarchical level, but the spins of the 3rd and higher 
level of hierarchy remain uncorrelated, and so on. Thus, 
approaching zero temperature is attended by a hierar- 
chy of phase transition, along which larger and larger 
clusters of spins become correlated. Note that the ap- 
plication of the spin-glass models to networks (see, for 
example, [17j ) suggests the opposite behavior, i.e. ap- 
proaching zero temperature is attended by extension of 
correlations from larger to smaller clusters of the net- 
work nodes. In this respect, the hierarchical networks 
may constitute an alternative approach to modelling of 
natural networks. 



motifs as in a special class of natural networks of essential 
biological importance. 

Let us point out two important features of hierarchi- 
cal networks constructed in our paper. First of all, any 
sub-graph belonging to a particular hierarchy level is an 
independent graph because the formation of clusters of 
links on each hierarchy level is entirely uncorrelated. Sec- 
ondly, the sub-graphs, associated with different hierarchy 
levels of the network, can be different, so the network as a 
whole can be essentially nonuniform. In nature, the ran- 
dom graphs of such a hierarchical genesis can be encoun- 
tered among the networks whose origins are associated 
with random events with low correlation, occurring with 
short evolutionary memory. The construction of such 
networks in some sense is very "simple", "rough" and 
"stable" because it does not demand a "fine tuning" of 
parameters to demonstrate the desirable properties (for 
example, the distribution of motifs). In particular, the 
networks of hierarchical genesis may by interesting as re- 
gards prebiology or the earliest biology. 

We believe that our result could shed the light on the 
relation between the distribution of motifs and the struc- 
ture of the adjacency matrix of a hierarchical network. 
However to make this relation more profound the "in- 
verse" problem should be considered as well. Namely, 
it would be desirable to check if the stable distribution 
of motifs is uniquely related to any kind of hierarchical 
organization of the network. The result of our work con- 
cerning the critical behavior of motifs' distribution on the 
noise intensity may be considered as a step towards this 
direction. This result demonstrates that the motifs's dis- 
tribution for hierarchical network has a "basin of stabil- 
ity" below some critical value / cr of random perturbation 
of the hierarchical network. 

We are guided by a general conjecture that the mo- 
tifs distributions corresponding to four supcrfamilies of 
U. Alon et al could signalize the existence of islands of 
stability (attractors) in a sea of possible motifs' distribu- 
tions. Whether this conjecture is true or not we hope to 
see in the close future. 



IV. CONCLUSION 



We have shown that the distribution of motifs in 
random hierarchical networks defined by nonsymmetric 
random block-hierarchical adjacency matrices coincides 
with the distribution of motifs in the second superfamily 
in the classification of U. Alon et al of networks [1JJ to 
which the class of neuron networks belongs. We would 
like to emphasize that apparently, this is the first ex- 
ample of "handmade" networks with the distribution of 
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Appendix A: Phase transition in hierarchical spin 
system in an external field 



In order to demonstrate that the existence of the phase 
transition is very natural in a hierarchical system, we 
consider the toy model of Ising spin system with symmet- 
ric block-hierarchical matrix of coupling constants. This 
model is known as "Dyson hierarchical model" J14l . K15j . 
We outline the standard derivation of the mean-field so- 
lution of a spin system and discuss briefly the obtained 
results for our specific hierarchy of coupling constants. 

Define the partition function Z of one-dimensional 
Ising spin chain with an arbitrary matrix U of coupling 
constants ua 



Z= £ e 

{si,...,sjv} 



i J2 U ijSi S j+Yl h i 



(Al) 



The spins s% (i = 1, ..., JV) take the values ±1 and hi 
is the external field acting on the spin Si. Using the 
Hubbard-Stratonovich transform 



±,J2 u ij s i s j t:~ n / 2 



VdetC/ 



N 



~\dxi 



-AT ^2wijXiXj J r^ J SiXi 



(A2) 



where wu is the element ij of the matrix W = U and 



substituting (|A2|) into (|A1|) . we can integrate over all spin 
configurations. The partition function Z reads now 



)-V 



\/ir N dct U 



A 



~[dx t 



-TF(xi,...xi 



(A3) 



where 



F(xi,...Xn) = iy]wjjXjXj 



1 N 

i=l 



In cosh(a;i 



(A4) 
We evaluate (|A3[) - (|A4[) in the saddle-point (mean-field) 
approximation. The partition function of the system is 

Z = cxp{-TF(x { °\ ...,x { ° } )} with F{x ( ° ] J0) 



by 

OF 



J N 



) given 
with Xi being the solutions of the equations 

0. Setting the relation between vectors 



{xl 



(0) 



x (0) 

[/-l x W ] we g e ^ : 



.((■») 

■N 



} and y 



(0) 



(0) 



{i/r,-,i/^}: y 



(0)i 



,{0) 



4y! 0) = itanh^ 



u ijVi + hi 



(A5) 



For small arguments of tanh(...) in the r.h.s. of (|A5P one 
can linearize (IA5I) and rewrite it as 



(A6) 



A' 



YJVii - 4T %) Vj + h i = o 
J'=l 



Suppose now that the matrix of coupling constants 
U = {uij} has the block-hierarchical structure identical 
to the structure of the Parisi matrix A shown in Fig|3] 
For simplicity we consider the matrix elements of A to 



be nonrandom with a 



(») 



iW 



2 -(«+i)7 ) an d t he ex- 
ternal field uniform hi = h for all i E [1, JV]. 

The solution of (|A6|) can be found using the methods 
of mathematical analysis on the field of p-adic numbers 
Q p . This technique is based on parametrization of the 
matrix elements \ug\ by the pairs of rational numbers 
{zi, Zj} (see [III Qi|) by such a way that the p-adic norm 
\zi — Zj\ p = pT( Zi - z i) induces block-hierarchical structure 
of the matrix U = {utj } and therefore the coupling con- 
stants Uij can be represented by an appropriate (real- 
valued) function Uij = p(\zi — Zj\ p ) — a^/ z . z .y In our 
case, p(\zi — Zj\ p ) is chosen in the form: 

fo, 7(*i,*j)<0 

p(\zi-Zj\ p )= < r iw ^ (A7) 

I 2 -(«+i)t(^,^) ! 1 (z i ,z j )>0 

Now, replacing the vector y(°) = {y[ ',..., y]^'} by the 
function f(zi) = y\ , we can read (|A6|) as follows: 

A 

£>(|* - ^-|p)/(^) - 4T/(z,) + h = (A8) 



Continuous analog of Eq. (|A8|) in the thermodynamic 
limit N — > oo is 



p(|z - z'| p )(^(2) d p z' - 47>(z) + /i = 0, (A9) 



Qp 



where z £ Q p and d p z is the Haar measure on 
Q p . Thus Eq. (|A8[) is understood as a discrete form 



of the p-adic equation (|A9[) induced by the relations 

J, z , <1 <yj(z) dj,z = f(zi) in the coset space Q p /Z p . 



The solution of Eq. (|A9[) is easily found using the p- 
adic Fourier transformation (see, for example, [16(). In 
terms of the Fourier transforms !p(k) and p(|fc| p ) of the 
functions (p(z) and p(\z — z'\ p ), Eq. (|A9p looks as 



p(\k\ p )ip(k) - ATip(k) + h = 
giving the solution 



AT-p{\k\ p ) 



(A10) 



(All) 



The function tp(k) has the poles at the set of critical 
temperatures T cr , determined by the equation 

4T cr - p{\k\p) = 

The values of T cr can be easily determined in the closed 
form by knowing that 



P(\k\ P ) 




\k\ P >l 
\kL<l 



(A12) 



7 



- i-p~ 



is the p-adic r -function and 



where r p (— a) - 

Since the Fourier transform (|A12[) of the function 
p(\z — z'\ p ) given by (|A7|) possess discrete values |fc|a = 



2 Q7 , 7 = 0, 1,2,..., the system of hierarchically interact- 
ing Ising spins has a hierarchy of critical temperatures 

T™ = p(2" Q7 )/4 (A13) 



It should be noted that in fact the hierarchy of critical 
temperatures T c ( r 7) defined by (|A12j) - (|A13|I exist only 
for a > 0. As soon as a — ► the intervals between T C r 
tend to zero and the hierarchical model becomes similar 
to ordinary mean-field ferromagnetic system without the 
hierarchy of interactions. 
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